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Abstract 

We consider a portfolio optimization problem in a defaultable market 
with finitely-many economical regimes, where the investor can dynami- 
cally allocate her wealth among a defaultable bond, a stock, and a money 
market account. The market coefficients are assumed to depend on the 
market regime in place, which is modeled by a finite state continuous time 
Markov process. We rigorously deduce the dynamics of the defaultable 
bond price process in terms of a Markov modulated stochastic differential 
equation. Then, by separating the utility maximization problem into the 
pre-default and post-default scenarios, we deduce two coupled Hamilton- 
Jacobi-Bellman equations for the post and pre-default optimal value func- 
tions and show a novel verification theorem for their solutions. We obtain 
explicit optimal investment strategies and value functions for an investor 
with logarithmic utility. We finish with an economic analysis in the case of 
a market with two regimes and homogenous transition rates, and show the 
impact of the default intensities and loss rates on the optimal strategies 
and value functions. 
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1 Introduction 



Continuous time portfolio optimization problems are among the most widely 
studied problems in the field of mathematical finance. Since the seminal work of 



Merton (1969), who explored stochastic optimal control techniques to provide 
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a closed form solution to the problem, a large volume of research has been done 
to extend Merton's paradigm to other frameworks and portfolio optimization 
problems (see, e.g. Karatzas et al. (1996), Karatzas and Shreve (1998), and 



Fleming and Pang ( 2004 )). Most of the models proposed in the literature rely on 



the assumption that the uncertainty in the asset price dynamics is governed by a 
continuous process, which is typically chosen to be a Brownian motion. In recent 
years, there has been an increasing interest in the use of regime switching model 
to capture the macro-economic regimes affecting the behavior of the market. 
More specifically, the price of the security evolves with a different dynamics, 
typically identified by the drift and the diffusion coefficient associated to the 
macro-economic regime in place. Although regime switching models arguably 
are able to incorporate a realistic description of market behavior, they pose 
challenges in the context of pricing because they lead to an incomplete market, 
as the regime uncertainty cannot be hedged away. 



In the context of option pricing, Guo (2001) studies a two-regime switch- 
ing model, where regimes represent the amount of information available to the 
market. [BufRngton and Elliott] and 'BufEn gton and Elliott] ( l2002-b I 

price European and American options under regime switching models, while 



Elliott et al. ( 2005 1 address the specification of an appropriate pricing martin- 



gale measure. Guo and Zhang (2004) provide an explicit optimal stopping rule 
when pricing perpetual American put options, while Guo et al. (2005j) consider 
the relation between regime shifts and investment decision in the context of 
real options. Graziano and RogersI (2006) give a methodology to price barrier 



options with regime switching dividend process, while |Gapeev and Jeanblanc| 
(2010) obtain closed form expressions for European claims, assuming geometric 
Brownian motion dynamics with regime switching drifts. 

Utility maximization problems under regime switching have been inves- 
tigated in Sotomayor and Cadenillas (2009), who consider the infinite hori- 
zon problem of maximizing the expected utility from consumption and ter- 
minal wealth in a market consisting of multiple stocks and a money market 
account, where both short rate and stock diffusion parameters evolve according 
to Markov-Chain modulated dynamics. Similarly, Zariphopoulou ( 1992 ) con- 
siders an infinite horizon investment-consumption model where the agent can 
consume and distribute her wealth across a risk-free bond and a stock. |Nagai| 



and Runngaldier (2008) consider a finite horizon portfolio optimization prob- 



lem for a risk averse investor with power utility, assuming that the coefficients 
of the risky assets in the economy are nonlinearly dependent on the Markov- 
chain modulated economic factors. Korn and Kraft (2001 ) relax the assumption 
of constant interest rate and derive expressions for the optimal percentage of 
wealth invested in the money market account and stock, under the assumption 
of a diffusive short rate process with deterministic drift and constant volatility. 

Most of the research done on continuous time portfolio optimization has 
concentrated on markets consisting of a risk- free asset, and of securities which 
only bear market risk. These models do not take into account securities car- 
rying default risk, such as corporate bonds, even though the latter represent 
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a significant portion of tlie market, comparable to tlie total capitalization of 
all publicly traded companies in the United States. In recent years, portfolio 
optimization problems have started to incorporate defaultable securities, but 
assuming that the risky factors are modeled by continuous processes and more 



specifically by Brownian Ito processes. Bielecki and Jang (2006) derive optimal 



finite horizon investment strategies for an investor with CRRA utility function, 
who optimally allocates her wealth among a defaultable bond, risk-free account, 
and stock, assuming constant interest rate, drift, volatility, and default inten- 



sity. Bo et al. (20101 consider an infinite horizon portfolio optimization problem. 



where an investor with logarithmic utility can choose a consumption rate, and 
invest her wealth across a defaultable perpetual bond, a stock, and a money 
market. They assume that both the historical intensity and the default pre- 



mium process depend on a common Brownian factor. Unlike Bielecki and Jang 



(2006), where the dynamics of the defaultable bond price process was derived 



from the arbitrage-free bond prices. Bo et al. (20101 postulates the dynamics of 



the defaultable bond prices partially based on heuristic arguments. |Lakner and| 
Liang ( 2008 ) analyze the optimal investment strategy in a market consisting of 



a defaultable (corporate) bond and a money market account under a continuous 
time model, where bond prices can jump, and employ duality theory to obtain 



the optimal strategy. Callegaro, Jeanblanc, and Runggaldier (2010) consider a 



market model consisting of several defaultable assets, which evolve according to 
discrete dynamics depending on partially observed exogenous factor processes. 



Jiao and Pham (2010) combine duality and dynamic programming to optimize 



the utility of an investor with CRRA utility function, in a market consisting of 
a riskless bond and a stock subject to counterparty risk. Bielecki et al. (2008) 



develop a variational inequality approach to pricing and hedging of a defaultable 
game option under a Markov modulated default intensity framework. 

In this paper, we consider for the first time finite horizon dynamic portfolio 
optimization problems in defaultable markets with regime switching dynamics. 
We provide a general framework and explicit results on optimal value functions 
and investment strategies in a market consisting of a money market, a stock. 



and a defaultable bond. Similarly to Sotomayor and Cadenillas (2009), we allow 



the short rate and the drift and volatility of the risky stock to be all regime 
dependent. For the defaultable bond, we follow the reduced form approach to 
credit risk, where the global market information, including default, is modeled 
by the progressive enlargement of a reference filtration representing the default- 
free information, and the default time is a totally inaccessible stopping time with 
respect to the enlarged filtration, but not with respect to the reference filtration. 
We also make the default intensities and loss given default rates to be all regime 
dependent. The use of regime switching models for pricing defaultable bonds 
has proven to be very flexible when fitting the empirical credit spreads curve of 
corporate bonds as illustrated in, e.g., Jarrow et al. (1997), where the underlying 



Markov chain models credit ratings. 

Our main contributions are discussed next. First, we rigorously derive the 
dynamics of the defaultable bond under the historical measure from the price 
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process, defined as a risk neutral expectation. Secondly, after separating the 
utility maximization problem into a pre-default and post-default dynamic opti- 
mization problem, we give and prove verification theorems for both subproblems. 
We show that the regime dependent pre-default optimal value function and bond 
investment strategy may be obtained as the solution of a coupled system of non- 
linear partial differential equations (satisfied by the pre-default value function) 
and nonlinear equations (satisfied by the bond investment strategy), each corre- 
sponding to a different regime. Moreover, we obtain the interesting feature that 
the pre-default optimal value function and bond investment strategy depend on 
the corresponding regime dependent post-default value function. Thirdly, we 
demonstrate our framework on the concrete case of an investor with logarith- 
mic utility, and, show that both the optimal pre-default and post-default value 
functions amount to solving a system of ordinary linear differential equations, 
while the optimal bond strategy may be recovered as the unique solution of a 
decoupled system of equations, one for each regime. Under a two-regime mar- 
ket with homogenous transition rates, we are able to obtain explicit formulas, 
and illustrate the impact of default risk on the bond investment strategy and 
optimal value functions. 

The rest of the paper is organized as follows. Section [2] introduces the mar- 
ket model. Section |3] derives the dynamics of the defaultable bond under the 
historical measure, starting from the risk-neutral bond price process. Section [4] 
formulates the dynamic optimization problem. Section [s] gives and proves the 
two verification theorems associated to the post-default and pre-default case. 
Section |6] specializes the theorems given earlier to the case of an investor with 
logarithmic utility. Section [7] summarizes the main conclusions of the paper. 
The proofs of the main theorems and necessary lemmas are deferred to the 
appendix. 

2 The Model 

Assume (ri,5,G,P) is a complete probability space, where P is the real world 
probability measure (also called historical probability), G = {Gt) is an enlarged 
filtration given by = TtW Tit (the filtrations J^t and Ht will be introduced 
later). Let {Wt} be a standard Brownian motion on {^,Q ,¥,V), where F := 
{J^t)t is a suitable filtration satisfying the usual hypotheses of completeness and 
right continuity. We also assume that the states of the economy are modeled 
by a continuous-time Markov process {Xt} defined on {n,Q ,¥,P) with a finite 
state space {xi,X2, ■ ■ ■ ,xn}- Without loss of generality, we can identify the 
state space of {Xt} to be a finite set of unit vectors {ei, 62, . . . , ejv}, where 
Bi = (0, 1, ...0) S and ' denotes the transpose. We also assume that 
{Xt} and {Wt} are independent. Define A{t) to be the Markov chain transition 
matrix [aij{t)]ij=i^2,...,N (also referred to as the infinitesimal generator). The 
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following semi-martingale representation is well-known (of. [Elliott et al. ( 1994 )): 

Xt^X^+ f A'{s)X,ds + Ivf{t), (1) 
Jo 

where M'^{t) = (Aff (t), . . . , Mj^{t)y is a IR^-valued martingale process under 
P, and A{t) — [a^ (t)]i j=i....^Ar is the so-called generator of the Markov process. 
Specifically, denoting Pi,j{t, s) :— V{Xs — j\Xt = i), for s > t, and Si j — li=j, 
we have that 

a,,(t) = lnn ^ ; 

cf. 



Bielecki and Rutkowski 



(2001). In particular, aisit) :~ 

We consider a frictionless financial market consisting of three instruments: 
a risk-free bank account, a defaultable bond, and a stock. The dynamics of 
each of the following instruments will depend on the underlying states of the 
economy as follows: 

Risk- free bank account. The instantaneous market interest rate at time t is 
rt :— r{t,Xt) :~ {r,Xt), where (•, •) denotes the standard inner product in 
and r = (ri, 7'2, . . . , r^)' are positive constants. This means that, depending on 
the state of the economy, the interest rate rt will be different; i.e., if = then 
Tt = Ti. The dynamics of the price process {Bt\ which describes the risk- free 
bank account is given by 

dBt = nBtdt. (2) 

Stock price. We assume that the stock appreciation rate {nt} and the volatility 
{(Tt} of the stock St also depend on the economy regime in place Xt in the 
following way: 

^lt := Xt) := (m, Xt) , at := Xt) := {a, Xt) (3) 

where = (/zi, /Z2, . . . , /iTv)' and a = {01,02, . . . , <Jn)' are constants denoting, 
respectively, the values of drift and volatility which can be taken depending on 
the different economic regimes. Hence, we assume that 

dSt = ^itStdt H- atStdWt, Sa = s. (4) 

Risky Bond price. Unlike the previous two securities, where we have written 
directly the dynamics under the historical measure, here we need to infer the 
historical dynamics (i.e. dynamics under the actual probability measure P) 
from the bond price process, which is originally defined under a suitably chosen 
risk- neutral pricing measure Q. Before defining the bond price, we need to 
introduce a default process. Let r be a nonnegative random variable, defined 
on C/, P), representing the default time of the counterparty selling the bond. 
Let Jit — a{H{u) : u < t) he the filtration generated by the default process 
H(t) :— lr<t, after completion and regularization on the right, and also let 
G := {Qt)t be the filtration Qt := FtW T-Lf 
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We use the canonical construction of the default time r in terms of a given 
hazard process {ht}t>o, which will also be assumed to be driven by the Markov 
process X. Specifically, throughout the paper we assume that ht :— {h,Xt), 
where h := {hi,h2, ■ ■ ■ jh^)' are positive constants. For future reference, we 
now give the details of the construction of the random time r. We assume the 
existence of an exponential random variable x defined on the probability space 
{^1,J-,P), independent of the process {Xt)t- We define r by setting 



: inf {t e 1R+ : /" hudu > 
Jo 



x}. 



(5) 



It can be proven that {ht)t is the {¥,G)-hazard rate of r (see Bielecki and| 
Rutkowski (20011, Section 6.5 for details). That is, {ht)t is such that 



-.^Hit)- f {l-H{u-))Kdu 
Jo 



(6) 



is a G-niartingale under P, where H{u ) = lims^„ H(s) — 1t<u- 

An important consequence of the previous construction is the following 
property. Let us fix t > and Too = Vs>o-^s- ^^'^ ™y " ^ have 
P(t < m| J'oo) = l-e- -fo Therefore, for any u<t, 

P(r < u\J^t) = [P(r < u\Too)\H = 1 " e-/o''*»'i- ^ p(r < u\Too)- 
Plugging u = t inside the above expression, we obtain 

V{T<t\Tt)=V{T<t\J^oc)- (7) 

It was proven in Bremaud and Yor] ( 1978 ) that Eq. ^ is equivalent to saying 
that any F-square integrable martingale is also a G-square integrable martingale. 
The latter property is also referred to as the H hypothesis, and we will make 
use of this property later on. 

The final ingredient in the bond pricing formula is the recovery process {zt)t, 
an F-adapted right-continuous with left-limits process to be fully specified below. 
Then, the time-t price of the risky bond with maturity T is given by 



^ e-~^>'''^''zudH{t 



^ Jt 



'{l-H{T)) 



(8) 



where Q is the equivalent risk-neutral measure used in pricing. Furthermore, 
we adopt a pricing measure Q such that, under Q, W is still a standard Wiener 
process and X is a continuous-time Markov process (independent of W) with 
possibly different generator A^{t) := [af^ (i)]i j=i,2,...,Ar- 

The existence of the measure Q in the previous paragraph follows from the 
theory of change of measures for denumerable Markov processes (see, e.g.. Sec- 



tion 11.2 in Bielecki and Rutkowski (2001)). Concretely, for i ^ j and some 



bounded measurable functions k^. 



— (— 1, oo), define 



4j{t) :=a,,,(t)(l + «;,,,(t)). 



(9) 
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and for i ^ j, define 



JV 



We also fix Ki^i{t) = for i = 1, . . . , iV. Now, consider the processes 

ft 



Ml'' := Hi'' 







(10) 



where 



Hl:=l{x,^,^}: and ff^'^ 



0<u<t 



.}l{X„=e,}, (*7^j)- (11) 



The process {Ml'-')t is known to be a F-martingale for any i j (see Lemma 
11.2.3 in Bielecki and Rutkowski (2001 1) and, since the i/-hypothesis holds 
in our default framework, they are also G-martingales. Then, by virtue of 
Proposition 11.2.3 in Bielecki and Rutkowski| (2001 ), the probability measure Q 
on G = {Gt)t with Radon-Nikodyn density {r]t} given by 



N 

1+ / E Vu-K^,J(u)dMl^^, 



(12) 



is such that X is a Markov process under Q with generator [afj{t)]ij=i,2,...,N- 
Without loss of generality, Q can be taken such that W is still a Wiener process 
independent of X under Q. 



(13) 



where is a E, -valued martingale under Q. In particular, note that 



M^{t) = M^{t)+ / {A{sy - A^{sy)X,ds. 



(14) 



We emphasize that the distribution of the hazard rate process ht = {h, Xt) under 
the risk-neutral measure is different from that under the historical measure. 
Therefore, our framework allows modeling the default risk premium, defined as 
the ratio between risk-neutral and historical intensity, through the change of 
measure of the underlying Markov chain. 
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3 Defaultable bond price dynamics 



We now proceed to obtain the bond price dynamics under both the risk-neutral 
and historical probability measures. Eq. ([8| may be rewritten as 



p{t,T) = 



(15) 



which follows from Eq. ([6|, along with application of the following classical 
identity 



Gt 



where t < s and F is a J^„-measurable random variable (see Bielecki and 



Rutkowski (20011, Corollary 5.1.1, for its proof). 



We assume the recovery-of-market value assumption, i.e. zt := {l—Lt)p{t^,T), 
where Lf is F-predictable. As with the other factors in our model, we shall 
assume that Lt is of the form Lf :— {L,Xt) for some constant vector L := 
(Li, . . . , Ln)' G K^. Under the recovery-of-market value assumption, it follows 
using a result in Duffie and Singleton] ( 1999 1, Theorem 1, that 



p{t,T) = lr>t 



The following result gives the dynamics of the defaultable bond price process 
under the risk-neutral measure Q. 

Theorem 3.1. Suppose that, for any i ^ j , the function afj defined in ^ is 
continuously difjerentiable in (0, T) such that 



< inf |af,-(s)| < sup jaf , (s)| < oo & sup 

s(^lPi:T] 'J se[0,T] se(0,T) 



ds 



< oo. (17) 



Then, the pre-default dynamics of the bond price p{t, T) under the risk-neutral 
measure Q is given by 



dp{t, T) = p{t , T) i [rt + ht [Lt -l)\dt+ ^^^^^ ^^^^ ~ d£,Y 



(18) 



where {M^{t))t is the N -dimensional (¥,Q) -martingale defined in (IS), {^'f)t 
is a {G,Q) -martingale, and '>p(t) :— {ipi{t), . . . ,'0Ar(i))' is given by 



Xt = Ci 



(19) 
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The proof is reported in Appendix]^ We also have the following dynamics 
under the historical probability measure. 



Proposition 3.2. Under the assumptions of Theorem \3.1\ the pre-default dy- 
namics of the bond price p{t, T) under the historical measure P is given by 



dp{t, T) = p{t-,T) I [rt + ht{Lt - 1) + D{t)] dt + ^V, ;, ^ - d^f 



'(20) 

where {M^{ty)t is the N -dimensional (F,P) martingale defined in (M), (^f)t is 
the {G,¥) -martingale defined in and D{t) {(Di(t), . . . , Dj^(t)Y , Xf) with 

4 Optimal Portfolio Problem 

We consider an investor who wants to maximize her wealth at time i? < T by 
dynamically allocating her financial wealth into (1) a risk-free bank account, (2) 
a risky asset, and (3) defaultable bond. The investor does not have intermediate 
consumption nor capital income to support her purchase of financial assets. 

Let us denote by the number of shares of the risk-free bank account that 
the investor buys {lyj^ > 0) or sells {lyf < 0) at time t. Similarly, ly^ and i^^ 
denote the investor's portfolio positions on the stock and risky bond at time i, 
respectively. The process {i/f ,1/^ ji/f) is called a portfolio process. We denote 
Vt{i') the wealth of the portfolio process = {i^^ , ^ v^) at time i, i.e. 

Vtiv) = lyfBt + vfSt + v^p{t.T). 

As usual, we require the processes , z/f , and to be F-predictable. We also 
assume the following self-financing condition: 

dVt = vfdBt + lyf'dSt + i^fdpit, T). 

Given an initial state configuration (x, z, u) G E := {ei, 62, ... , e^} x {0, 1} x 
(0,00), we define the class of admissible strategies A :— A{v,i^z) to be a set 
of (self-financing) portfolio processes v such that Vt{v) > for alH > when 
Xq — X, Hq = z, and Vq ~ v. 

Let TT := (tt^, 7rf , TT^) be defined as 

if Vt^{iy) > 0, while irf = 7rf ^ vrf = 0, when Vt-{iy) = 0. The vector tt := 
{nf , 7rf , Trf ), called a trading strategy, represents the corresponding fractions of 
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wealth invested in each asset at time t. Note that if tt is admissible, then the 
dynamics of the resulting wealth process can be written as 



dBt 
Bt 



St 



dp{t,T) 
p{t-,T) 



under the convention that 0/0 = 0. This convention is needed to deal with the 
case when default has occurred (t > r), so that p{t~ ,T)=0 and we fix 7rf = 0. 
Using the dynamics derived in Proposition 3.2 and that tt^ + tt 



TT^ = 1, we 



have the following dynamics of the wealth process 



^t- 



{rt + 4{[it - rt) + T^^{ht{Lt - 1) + D{t)W dt 



{m,dM^{t)) 



(23) 



under the actual probability . 



4.1 The utility maximization problem 

For an initial value {x,z,v) £ E and an admissible strategy tt — (tt^, tt^, tt^) S 
A{x, z, v), let us define the objective functional to be 



Jr{x,z,v;'k) := 



Xq = x,Ho = z,Vo ^ V 



(24) 



i.e. we are assuming that the investor starts with v dollars (its initial wealth), 
that the initial default state is z (z = means that no default has occurred 
yet), and the initial value for the underlying state of the economy is x. The 
constraint Vq = v is also called the budget constraint. As usual, we assume that 
the utility function U : [0, oo) — !■ M U {oo} is strictly increasing and concave. 

Our goal is to maximize the objective functional J{x, z,v;tt) for a suitable 
class of admissible strategies ttj := {irf ,nf ^nf). Furthermore, we shall focus 
on feedback or Markov strategies of the form 

TTt^in^ [t,Vt-,H{t-),Til {t,Vt-,H{t-),TTP {t,Vt-,H{t-)), 



for some functions 7rf , 7rf , 7rf : [0, oo)x [0, oo)x{0, 1} — > M such that 7r^(t, w, z)+ 
TTf{t,v,z)+Tif{t,v,z) = 1. 

As usual, we consider instead the following dynamical optimization problem: 



ip^{t,v,i,z) := sup E"" 



Vt^v,Xt^e,,Hit) = z 



(25) 
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for each (w, i, z) G (0, oo) x {1, 2, . . . , N} x {0, 1}, where 



{r, + Trf (^i, - rs) + Trf (1 - H{s-))[hs{Ls - 1) + 



+ 7rfa,dVF, + 7rf(l-H(.s-)) 



(V.(.),dAfP(,s)) 



VV =1). 



(26) 



The class of processes At{v,i,z) is defined as follows: 



Definition 4.1. Throughout, At{v,i, z) denotes a suitable class of ¥ -predictable 
locally bounded feedback trading strategies 

TT, (^f,7rf) (vrg^ Js, ^/j*^", (s")), _ (s, i7(s-))), se[t,R], 



such that (26) admits a unique strong solution {V7^'*'"}se[t,_R] and Vg'*"'"" > 
for any s £ [t, R] when Xt = and H{t) = z. Throughout this paper, a trading 
strategy satisfying these conditions is simply said to be t-admissible (with respect 
to the initial conditions Vt = v, Xt = e^, and Ht = z). 



Remark 4.1. As it will be discussed below (see Eqs. (31), (32), and {6(^ ), the 

jump AK, Vs ~ Vg- of the process (26) at time s is given by 



N 



i=l j^i 



AHl'^^TT^AHis) . (27) 



Since for {Vs}s>t to be strictly positive, it is necessary and sufficient that AVg > 
—Vg- for any s > t a.s. (cf. ^Jacod and Shiryaev 2003, Theorem 4-61)), we 
conclude that in order for vr^ to be admissible, it is necessary that. 



M, 



^i(s) 



< 7r,f (s, V, z) < 1, 



(28) 



for any s,v > 0, z & {0,1}, and i = 1,...,N, where we set Mi :— — oo if 
ipi{s) > ipj{s) for all j ^i. 



5 Verification Theorems 

As it is usually the case, we start by deriving the HJB formulation of the value 



function ( 25 ) via heuristic arguments. We then verify that the solution of the 
proposed HJB equation (when it exists and satisfies other regularity conditions) 
is indeed optimal (the so-called verification theorem). Let us assume for now 
that ip^{t, V, i, z) is in t and in v for each i and z. Then, using Ito's rule 
along the lines of Appendix |Bj we have that 

^«(^,F,^C^,i/(^)) =v^^(r,y/,a,i7,)+ / C^''{s,V: ,Cs,Hs)ds + Mt-Mr, 
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where {Ct)t is the Markov process defined in (62 1, £ is the infinitesimal generator 
of (t, Vt,Ct,Ht) given in Eq. (69), and {Mt)t is the martingale given by Eq. (70). 
Next, if r < t < i?, by virtue of the dynamic programming principle, we expect 
that 

(^^(r, y/, Cr,Hr) = maxE [</7«(i, V^^ , Ct, HmSr] ■ (29) 



Therefore, we obtain E 



Gr 



< 0, with the inequality 



becoming an equality ii n — tt, where tt denotes the optimum. Now, evaluating 
the derivative with respect to t, aA, t — r, we deduce the following HJB equation: 



max Cip^' {r, V, i, z) = 0, 



(30) 



with boundary condition (p^{T, v, i, z) = U (v). 



In order to further specify ( [30| , let us first note that the dynamics ( 26 ) can 
be written in the form 



N 



dVs=a^Js + 'd^JW, + J2Pc _,,dMj{s)--f^ _df^, {t<s<R), (31) 



with coefficients 



Pi,j{t,v,z) = VT:[{t,v, z){l - z) 



li{t,v,z) = VTr[{t,v,z){l - z), 



a,(t, v,z)^v [r, + Trf (^, - n) + Trf (1 - z){h,{L, - 1) + A(<))] (32) 
'di{t,v,z) = Tif {t,v,z)aiV, 



where Di{t) is defined as in ( [21) ). Using the expression for the generator in 
Eq. (69), the notation (pi^z{t,v) := (p^{t,v,i, z), and the relationship tt^ = 
1 — tP~^ tt^ , ( 30 ) can be written as follows for each i = 1, . . . ,N: 







dt 



dv 



max \ TT, (^i - n)v-^ + (tt, ) 



+ (1 - z) max I TT[0i{t)v—^ + hi [ipi^i{t,v{l - 7rf )) - (pi^z{t,v)] 



1 + n 



where 



We can consider two separate cases 



(33) 



(34) 



ip^{t, V, i) — fi^oit, v) — ip^{t, V, i, 0), (pre-default case) 



(35) 
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and 

ip^{t,v,i) — ipi^i{t,v) — ip^{t,v,i,l), (post-default case). (36) 



Sect ion [O] give a verification theorem for the post-default case, while Section 
5. 2 1 gives a verification theorem for the pre-default case. 



5.1 Post-Default Case 

In the post-default case, we have that p{t,T) — 0, for each t < t < T. Conse- 
quently, nf — tor T < t < T and, since nf = 1 — 7rf — nf, we can take tt — 
as control. 

Below, rji :— denotes the sharpe ratio of the risky asset under the 

i*'' state of economy and Cq'^ denotes the class of functions zu : [0, R] x x 
{1, . . . , TV} ^ M+ such that 

ro(-,-,i) e C^'2((0,i?)xM+)nC([0,i?]xM+), TUy{s,v,i)>0, w^y{s,v,i) < 0, 

for each i — 1, . . . ,N. We have the following verification result, whose proof is 
reported in Appendix [C| 

Theorem 5.1. Suppose that there exist a function w e Cq'^ that solve the 
nonlinear Dirichlet problem 

2 2 / -A 

/ ■ \ V ^ ; ^ ; ^ j / ■ \ \ ^ /\// ■\ / -\\ r\ 

m[s, v,^)-^ 7 T +r^vw^{s, v,i)+y a,.j (s) (w(s, v,j)- w{s, v, i)) = 0, 

(37) 

for any s € (0, R) and i ~ 1, . . . , N , with terminal condition w{R, v, i) = U{v). 
We assume additionally that w satisfies 



(i) \wis,v,i)\<Dis) + Eis)v, (ii) 



w^is,v,i) 



Wyvis,v,i) 



<Gis)il + v), (38) 



for some locally bounded functions D,E,G : M+ — > M+. Then, the following 
statements hold true: 

(1) w{t,v,i) coincides with the optimal value function {t , V , i) = ip^{t,v,i,l) 



in (25), when At{v,i,l) is constrained to the class of t- admissible feed- 



back controls 7rf = tt^^ (s, Vg) such that TTi{-, ■) € C'([0, R] x M+) for each 



i = 1, . . . ,N and 



\vMs,v)\<G{s){l + v), (39) 
for a locally bounded function G. If the solution w is non-negative, then 



condition {39) is not needed. 



(2) The optimal feedback control {T^s}s£[t,R), denoted by , can be written 
as 7rf — TT^s (s, Vs) with 

Tii{s,v) = -. (40) 

a^ VW^^{S,V,l) 
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5.2 Pre-Default Case 



In the pre-default case {z = 0), we take and tt^ as our controls. We then 
have the following verification result: 



Theorem 5.2. Suppose that the conditions of The orem 5.1 are satisfied and, 
in particular, let w G Cq'^ be the solution of (SI). Assume that w G Cq'^ 

N , solve simultaneously the following system of 



and Pj = pi{s,v), i = 1, 
equations: 

9i{s)'Wv{s, V, i) - hii£^{s, v{l - Pi), i) 



1 \ Wv [ S,V 



- 1 



,J =0, 



(41) 



- / -N wl{s,v,i) _ 

wt{s,v,i) - — + r^vw^(s,v,i) 

2 w^v{s,v,i) 



+ \ Pi0.i{t)vwv{s, V, i) + hi [w(s, vil - Pi), i) ~ w{s, v, i) 

>j(s) 



w [ s,v I 1 + p. 



i^iis) 



- 1 



j - w{s,v,i) 



= 0, (42) 



for t < s < R, with terminal condition w(R,v,i) ~ U{v). We also assume 
that pi{s,v) satisfies (28) and [39) (uniformly in v and i) and w satisfies (38). 
Then, the following statements hold true: 

(1) w{t, V, i) coincides with the optimal value function (p^{t, v, i) — <p^(t, v, i, 0) 
in ^25\), whenAt {v,i,0) is constrained to the class of t-admissible feedback 



controls {nf ,Trf) = [tt^ {s,Vs-,H{s )),tt^ {s,Vs~,H{s ))) such that 

7Tf{;;z),7T[{-,-,z)eCi[0,R]xR+), 



for each i = 1, . . . , N , tt^ sati sfies (39) for a locally bounded function G, 
and TT^ satisfies (28) and (39) (uniformly in v,i,z). If the solution w is 
non-negative, then these bound conditions are not needed. 

(2) The optimal feedback controls are given by Trf :— irf, {s,Vs, H{s)) and 



{t,Vt,H{s)) with 



7rf{s,v,z) 
Tr[ {s,v,z) 



rj, w^{s,v,i) ?7, w^{s,v,i) _^ ^^^^ 



CTj vwi,v{s, V, i) 
Pt{s,v){l - z) 



CFi vw^^{s,v,i) 



(44) 



6 Application to Logarithmic Utility 

The objective of this section is to specialize the framework developed in the 
previous sections to concrete choices of utility functions. We focus on the loga- 
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rithmic utility function U{v) = log(z;). The framework, however, can be appUed 
to other concave utiUty functions, such as power or negative exponential utilities. 
For the sake of clarity and conciseness, the details of the numerical implementa- 
tion for other HARA functions is being deferred to the follow-up paper [Capponi,| 



Figueroa-Lopez, and Nisen (20111. 



6.1 Explicit Solutions 

Let us recall that the investor's horizon R is assumed to be less than the ma- 
turity T of the defaultable bond. Before proceeding, we state without proof a 
fundamental result from the theory of ordinary differential equations. 



Lemma 6.1 (Codd (1961)). Suppose that the nx n matrix F{t) and the nx 1 
vector b{t), are both continuous on an interval / € JR. Let to d I . Then, for 
every choice of the vector xq, we have that the system 

x'{t)^ F{t)x{t) + b{t), x{ta)=xo, 

has a unique vector-valued solution x{t) that is defined on the same interval I. 

We first give a lemma, which will be used later to derive the dynamics of 
the optimal value functions, and relations satisfied by the optimal investment 
strategies. 

Lemma 6.2. The system of equations 

for i = l,...,iV, admits a unique real solution Pi(s) in the interval (Mi,l), 
where Mi G [— oo, 0) is defined as in (28). Moreover, if for each i,j — 1, . . . , iV, 
Qij and a^j are continuous functions, then p{s,i) is a continuous function of s. 

The proof of this lemma is reported in Appendix [Dj The following is our main 
result in this section. 

Proposition 6.3. Assume that the afj 's and ai j 's are continuous in [0, T]. Let 

p{t) ~ [pi{t) , p2{t) , . . . ,PN{t)] be the unique continuous solution in [0,T] of the 
nonlinear system of equations UW. Then, the following statements hold true: 



(1) The optimal post-default value function is given by 

ip'^ {t, V, i) =\og{v) + K{t,i), 

where < t < R, and K{t) ^ [K{t,l), K{t,2), . . . K{t,N)] is the unique 
solution of the linear system of first order differential equations 

Kt [t, i)+r., + ^+Y. (*)-^(^' ■^') + (^)^(^' ^ 0, (46) 
with terminal conditions K{T, i) = 0, for i = 1, . . . , N . 
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(2) The optimal percentage of wealth invested in stock is given by tt (t) = 
[nf{t),ni{t), . . .,irf,{t)], where (t) = 0<t<R. 

(3) The optimal percentage of wealth invested in the defaultable bond is 7ff = 

(t), while the optimal pre- default value function is 

ip^it,v,i)=\og{v) + J{t,i), 

where J{t) — [J{t, 1), J{t, 2), . . . J{t, N)] is the unique solution of the linear 
system of first order differential equations 

2 

Jt{t, i) + '^+r^+ p^{t)e,{t) + h, (log(l - p,{t)) + K{t, ^) - J{t, ^)) + 



log 



1 



with terminal conditions J{R, i) — 0, for i 



+ J{t,j)-J{t,i) 



0, (47) 



1, 



The proof of this proposition is reported in Appendix |D] We notice that the 
only difference between the pre-default and post-default optimal value function 
lies in the time and regime component. Moreover, we obtain that the optimal 
proportion of wealth invested in stocks is constant in every economic regime, 
and does not depend on the time or on the current level of wealth. This is 



consistent with the findings in Sotomayor and Cadenillas ( 2009 1 and Bo et al. 



(2010). We also find that the optimal proportion of wealth allocated on the 
defaultable bond is time and regime dependent, but independent on the current 
level of wealth. Bo et al. (2010) find that the optimal allocation only depends 
on time through the default risk premium. 

In the case where the rate matrix A is homogenous (i.e. aij{t) = Oij) and 
the number of regimes is = 2, we can obtain closed form expressions for 
the optimal pre-default value function, post-default value function, and optimal 
bond fraction tt^ . Let us introduce the following notation: 



Q = n + 



2 ' 



Mt) 



Cl = hi — aij, C2 — /i2 



Mt) 

- 0-2.2, 



1, 



Mt) 



Mt) 



1, 



C+ = (C2 - Ci)^ + 401402,2, 



g,it) = n[{t)9,{t) + h, (log(l - ^P{t, i)) + K{t, i)) - o,,, log (l + ^f{t)Mt)) 

A,{t) = ~Ae,(t)Mt){K - 6,{t) + auMt)) + {0,{t) - a^,,Mt) + h^Mt) - 0^(t)^^(^))^ 
for i ~ 1,2. We have the following result: 

Corollary 6.4. Assume N = 2 and the rate matrix A to be homogenous. Then, 
the solution functions K{t, 1) and K{t, 2) of {^6) are given by 

1 

(ai,i + 02,2) 

+Q^[4-^,3MR-t)+a,., (e^-'^'i 



K{t,i) = 



Ca-.OM (1 - e(-i.i+«^>^)(«-*) + (ai,i + a2,2)(i? - i)) 

l + 03-i,3-« {R-t)\ 



l+a-2,2)(R-t) 
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The optimal feedback bond-fractions functions (t) and {t) are given by 



hi+aiiil>i{t) 



20i{t)ii,(t) 



ife,{t) = Q, 
, ife,{t)^Q. 



Finally, the solution functions J{t, 1) and J{t, 2) of (4-1) are given by 
J{t,i) 



1 i(<=l+g2) 

2c+ 



-^c+ cosh ( ) + (c2 - ci) sinh 



R 



(48) 



X e ''^T'^^'' (2{e'V^ - l)(/i, - q)(C3-. + 53-*(s)) 
+ ((3 - 2z)(c2 - ci) - - e^^((3 - 2i){c2 - ci) + yH^))(G + g^{s)) ]ds 



2(/ii — Ci) sinh ( — I x / e 



(-V^ + e^^[(3 - 2z)(c2 - ci) - VcT] - (3 - 2i)(c2 - ci))(C3-^ + 9^-^{s)) 
+ 2(e^VST _ i)(/^3_^ _ C3_,)(G + 9^[s))]ds . 



From the formulas given above, we can notice that for any given value of the 
historical intensity ft-i, the optimal bond strategy Trf (i) is independent of the 
value of the historical intensity h2, associated with the other regime. A sym- 
metric argument applies to Trf (t), thus showing that each strategy depend on 
the other regime only through the risk-neutral loss and default intensity asso- 
ciated to the other regime. This is not surprising, given that an investor would 
base his decision to buy or sell a defaultable bond on the market perception of 
default risk (i.e. based on the risk-neutral default intensity parameters) rather 
than on the number of defaults experienced by the corporation in the past. 



6.2 Economic Analysis 

The objective of this section is to measure the impact of the default parameters 
over the value functions and the optimal bond strategy via numerical analysis. 
We fix the interest rate, drift, and volatility regime parameters as indicated 
in Table [TJ We choose the transition rates of the chain to be the same under 
both probability measures. In order to evaluate ipi(t), we use the analytical 
expression for the probability density fi (R, x) of the time spent in state i by 
a two-regime continuous time Markov Chain, for a given time interval [0,i?], 
when the chain starts in state i at time 0. Such formulas have been provided 



in Kovchecova et al. (2010). Applying these formulas to our setting, for a given 
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Regime '1' 


Regime '2' 


r 


0.03 


0.03 




0.07 


0.02 


a 


0.2 


0.2 



Table 1: Parameters associated to the two regimes. 




< a; < i?, we obtain 

where S{x) denotes the Dirac deha function and Ij(z) is a modified Bessel 
function of first kind. When applying this formula to our case, we have 

Jo 

R. 

-{ri+hiLi){R-x)-(r2+h2L2)i 







We plot the behavior of 7rf (0) and (0) with respect to hi and /i2. It 
appears from Figure [T] that the strategies are not very sensitive to the starting 
regime when the ratio is close to one, with hi and /12 not too large. This is 
consistent with intuition, because in such a scenario, the default intensity is the 
same and small under both regimes, thus the default probability is also small 
and, consequently, the slightly different risk neutral loss rates in Table [l] do not 
affect much the investment choice. As the gap between /12 and hi increases, 
the strategies 7rf (0) and 7f|'(0) behave differently. This is because a larger 
risk-neutral intensity translates into a larger risk-neutral default probability, 
and given that the risk-neutral transition rates of the chain are not large, the 
starting regime matters. 

We next show the behavior of the optimal bond strategy over time. It 
appears from the left panels of Fig. [2] that the number of bond units sold de- 
crease as the investment horizon increases. It can be noticed that, the riskier 
the corporate bond, the smaller the number of units sold. This is negatively 
correlated with the regime conditioned bond prices, as it can be checked from 
Fig. [2] showing that riskier bonds have smaller prices, and shorter maturity 
bonds have larger prices. Moreover, since the second regime is riskier and the 
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Figure 1: Optimal bond strategy vrf (0) and '?t|'(0) versus the risk neutral hazard 
intensities hi and /i2- The upper surface represents 7rf (0), while the lower 
surface represents 7r|'(0). The loss parameters are given, respectively by Li = 
0.4, and L2 = 0.45. The transition rates are 012 = 0.7, and 021 = 0.1. The 
horizon R is set to 2 years. 



transition rate from the second to the first regime very small, the price decay 
over time of the bond is faster when the Markov chain starts in the second 
regime (V'2(i) < thus leading to a larger number of bond units sold when 

the chain starts in the second regime. 

We next inspect the behavior of the time components, K{t, 1) and K{t, 2), of 
the optimal post-default value function over time, for a fixed investment horizon 
R. We report the results in Figure |3j and in each plot we superimpose the time 
component of the optimal value function obtained in the Merton model. The 



latter is well known from the work of Merton (1969), and in our specific case 



obtained assuming that for each regime i, the market model consists of a stock 
and a money market account with parameters /ii, r^, and ai. It appears from 
the plots that both post-default value functions decrease with time. Moreover, 
they differ at times t far from the investment horizon due to possibility of regime 
shifts, and start converging to each other when the time to horizon is small. This 
is because the chain spends the largest fraction of its time in the starting regime 
due to the small transition rates (012 — 0.4 and 021 = 0.1), and consequently 
for short times to the horizon, the wealth process of our regime switching model 
approaches the wealth process in the Merton model. 

We finally evaluate the behavior of the time components, J(t, 1) and J{t, 2), 
of the optimal pre-default value function over time, assuming again a fixed 
investment horizon R. We report the results in Figure [4| where we vary /12, 
keeping hi fixed. We notice from the plots that the pre-default value function 
is decreasing with time, and very sensitive to the default risk level. 
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0.5 



1 1.5 2 

Horizon R 



2.5 3 



0.5 1 1.5 2 2.5 3 

Horizon R 




Figure 2: Optimal bond strategies and regime-conditioned bond prices at time 
zero versus time horizon. The top left panel represents the optimal bond strat- 
egy Trf (0). The top right panel represents the regime-conditioned bond price 
ipiiQ)- The bottom left panel represents the optimal bond strategy TrflO). The 
bottom right panel represents the regime-conditioned bond price V'2(0). The 
loss parameters are given, respectively, by Li = 0.4, and L2 = 0.45. The hazard 
intensity hi = 0.1. The transition rates are 012 = 0.4, and 021 = 0.1. 
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Figure 3: The left panel represents the optimal post default value function 
K{t, 1), and the optimal value function in the Merton model with parameters 
/Zi , ri and CTi . The right panel represents the optimal post-default value function 
K{t,2), and the optimal value function in the Merton model with parameters 
yU2, r2 and a2- The transition rates are ai2 = 0.4, and 021 =0.1. The horizon 
R is set to three years. 




0.5 1 1.5 2 2,5 3 0.5 1 1,5 2 2.5 3 



Figure 4: The left panel represents the optimal pre-default value function J(t, 1) 
for different levels of default risk associated to regime 2, and parameterized 
by the pair (/i2,-^2)- The right panel represents the optimal pre-default value 
function J{t,2) for the same levels of default risk (ft.2,i2) as in the left panel. 
The transition rates are ai2 — 0.4, and 021 =0.1. The horizon R is set to three 
years. 
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7 Conclusions 



We considered the continuous time portfolio optimization problem in a default- 
able market, consisting of a stock, defaultable bond, and money market ac- 
count. We assumed that the price dynamics of the assets are governed by a 
regime switching model. We derived the dynamics of the defaultable bond un- 
der the historical measure from the risk neutral price process. We have shown 
that the utility maximization problem may be separated into a pre-default and 
a post-default optimization subproblem, and proven verification theorems for 
both cases under the assumption that the solutions are monotonic and concave 
in the wealth variable v. The post-default verification theorem shows that the 
optimal value function is the solution of a nonlinear Dirichlet problem with 
terminal condition. The pre-default verification theorem shows that the opti- 
mal pre-default value function and the optimal bond investment strategy can 
be obtained as the solution of a coupled system of nonlinear partial differential 
equations with terminal condition (satisfied by the pre-default value function) 
and nonlinear equations (satisfied by the bond investment strategy). Each equa- 
tion is associated to a different regime, and the dependence of a regime i from 
another regime j comes through the Markov transition rates and the ratio be- 
tween the defaultable bond prices in regime j and regime i. Our results imply 
that the pre-default optimal value function and the bond investment strategy 
depend on the optimal post-default value function. 

We demonstrated our framework on the concrete case of an investor with log- 
arithmic utility, and shown that both the optimal pre-default and post-default 
value function can be obtained as the solution of a linear system of first order 
ordinary differential equations, while the optimal bond strategy can be uniquely 
recovered as the solution of a decoupled system of nonlinear equations, one for 
each regime. We have also performed an economic analysis on a two-regime 
market model with homogenous transition rates, and investigated the impact 
of default risk on the optimal strategy and value functions. Our analysis has 
shown that the optimal number of bond units sold in each regime decreases 
with the riskiness of the bond perceived by the market, and that the number 
of bond units sold is smaller for larger investment horizons. Although we have 
specialized our framework to the specific case of logarithmic utility, it is flexi- 
ble enough to accommodate any concave increasing utility function. Therefore, 
the results derived in the verification theorems can be used to derive explicit 
or numerical solutions for pre-default and post-default value functions, as well 
as optimal investment strategies, corresponding to a wide range of utilities of 
practical interest, such as the ones in the HARA family. 



A Risk-neutral and historical bond dynamics 

The following result states that the process introduced in ^ is also a Q- 
martingale. However, in order to indicate in the sequel when certain dynamics 
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are being taken under Q or under P, we will introduce a new notation . We 
should keeping in mind through the proof below that ^f. 



Lemma A.l. The process 

^f:=H{t)- l\l~H{u-))K 
Jo 

is also a {G,Q)— (loeal) martingale. 



du 



(49) 



Proof. By the definition of Q and the fact that 77 is a G-martingale (see the 
paragraph before (12)), it suffices to prove that rjt^^ is ^ G-martingale under P. 
From Ito's formula and the definition of 77 in ( 12 1, we have the process 

Jo Jo 



From (|6|, (12), and (49), ir]t^f)t can be written as 



0<s<tk,l=l 



where H^'^ := J2o<s<t ^{x _=k}^{x^=i}- Since the first two terms on the right- 
hand side of the previous equality are (local) martingales under P, it remains 
to show that the last term vanishes. But, given that AHg 7^ at s = t, 
the summation in the last term above will be provided that AXr = 0, a.s. 
In order to show this, let us recall that by definition X has no fixed-jump 
times; i.e. P{AXt 7^ 0) = for any fixed time t > 0. Also, using the def- 
inition of T given in Eq. (jij), r = inf{< > : h{Xs)ds > ^}, where ^ is 
an exponential random variable independent of X. Then, conditioning on X, 
¥{AXr^O) =E{E[l{AX,#o}|^s,s>0]}. Denoting < n < T2 < . . . the 
transition times of the Markov chain X, E [ifAx^^^o}] ^s, s > O] is given by 



E 



Xs,s> 



E 



(Xs)ds=i} 



Xs,s>0 



where the last equality follows from the independence of X and ^, and the fact 
that ^ is a continuous random variable. □ 



Lemma A. 2. Under the assumptions of Theorem\3.1\ the function 



V^,(i) = EQ 



Xt 



(50) 



is differentiable for any t € (0, T). 
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Proof. Define functions kij : [0,oo) — )■ (— l,C)o) such that 

j^\l = + f^ijit))' ioi i^j, and k;^, = 0. 

We also let Oij :— 1/{N—1) for i ^ j and di^i — —1, so that A :— [a^.j] ■ ^-^^ ^ 
is a valid generator of a homogeneous Markov process with transition times 
determined by a homogeneous Poisson process and an embedded Markov chain 
{Xi}i>i with transition probabilities pij := 1/{N — 1) for i ^ j. Now, let us 
define a probability measure Q with Radon-Nikodyn density (fit)t given by 



f ^ ~ 
1+ / Y] f]u-~^^,J{u)dMl;\ 
J {OA ,^^.^1 



(51) 



where Ml''' := H^"-' — C af^{u)Hl^du, and we us ed notation By virtue 

of Proposition 11.2.3 in [Bielecki and RutkowsE (20011, {Xt}t>o is a contin- 



uous Markov process with generator A under Q. Next, note that iplt) 
{ipi{t), . . . , V'Ar(i)) is such that 



But also, changing into the probability measure Q, we can write 



(52) 



- J^(,rs+h,L,)ds 



and, hence, we have the following representation for -0^: 



VT 



Recall that the solution of (51 ) can be written as 



Vt 



Let K{t) = [K,^j{t)],^j and f{t) := . . . ,rjv(t))' be defined by 



N 



K,,j{t) := log (1 + k,^,{t)) , f,{t) := r, + h,L, - ^ a%{t)k,^,{t). 



Then, we have 



exp ■ 



f{s)'X,ds ■ 



^ K-K{s)X, 



Xt = ej 
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where X' denotes the transpose of X. Next, using that {Xt)t is a homogeneous 
Markov process under Q, 



exp ■ 



T-t 



r{t + s)X,ds- X'^-K{t + s)X, 

s6(0,T-t]:AXa5^0 



where we used the notation (•) := E'^ ( -j Xq = e^). Furthermore, in terms of 
the transition times Ci < C2 < • • • of X, the embedded Markov chain {Xi} of 
X, and the number Mt of transitions by time t oi X^ can be written as 



C,=T-t 



where Co = 0. Using that {Mt)t is a Poisson process under Q and conditioning 
on A/^, we have 



m=0 



' 11,1=0 IcL 



' f{t+sYX„ds-J2'^^i X'„-iK(t+C,U(^„))X„ 



where ( = T — t, [/(i) < C/(2) < • • • < f/(m) are the ordered statistics of m i.i.d. 
uniform [0, 1] variables independent of X, U(^q) = 0, and J7(m+i) = 1- From the 
previous expression, we see that it suffices to show that 



e 



(53) 

is continuously differentiable in C g (0, T) for each to > 0, and that there exists 
a sequence {K„i}m>o such that 



(i) sup |$m(C)|<i^m, (ii) sup |$™(C)|<ifm, (iii) V e ^^K.^<(X). 

0<C<T 0<«T 



m— 



We now show that (i)-(iii) are satisfied provided that 



sup |fi(s)| < 00, sup \Kij{s)\ < 00, and sup \Kl j{s)\ < 00, 

sG[0,T] se[0,T] ' se[0,T] 



are satisfied. The latter conditions directly follow from (17 1. Denoting TmiC) 



the random function inside the expectation E^ in (|53|), one can check that 

„vaT maxi su 



sup |r™(c)|<e^ 

o<C<T 



Po<<:<T l''i(C)'l+™maXi.3 supo<^<y |ifi.3(C)| J(m ^ (^54^ 

for a constant A < oo. Also, TmiC) is continuously differentiable and 

f m 

m '\ 

YK-lK'{T - CUin))XnUin)\ , 
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where C/(„) := 1 — ?7(„). In particular, there exists a constant B < oo such that 

sup \r'„^iO\ < BmA^r (55) 

0<C<T 

By the formal definition of the derivative $J„(C) and the dominated convergence 
theorem, one can check that ( 54||55 1 will suffice for (i)-(iii). □ 

We are ready to give the proofs of the bond price dynamics: 



Proof of Theorem 3. 1 We first write the pre-default dynamics of the bond price 
under the risk neutral measure (i.e. on the event t > t). From Eq. (16), we 
have 



Define v{t, T) 



^-J^{r,+KL,)ds]^xi\ and H{t) = lr>t so that p{t,T) 



Htv{t,T). In terms of (19 1, note that v{t,T) — {ip{t),Xt), where ip{t) is given 
as in (19). In particular, Av{t,T) = (V'(t),AX(), since ip{t) is continuous 
in light of Lemma A. 2 Next, let us introduce the (F, Q)-martingale 0(t) := 



e-Io irs+hsL,)ds^jr^ ^ g^^^ ^jjg process b{t) := exp{/Q(r^ + /i,L,)ds}. Then, 

v{t,T) = b(t)4){t) has dynamics dv(t,T) = {n + htLt)v(t,T)dt + b(t)d(l>{t). This 
leads to 

dp{t,T) = H{t-)dv{t,T) + v{t-,T)dH{t) + Av{t,T)AH{t). (56) 

By virtue of the identity Av{t, T) — {ip{t), AXt) and similar arguments to those 
in the proof Lemma |A.1[ we have that 

J2 Avis,T)AH{s) ^ J2 (i'.,^XsAH{.s)) ^ {^j;r,AXr) =0. 



0<s<t 



Q<s<t 



Thus, (56) simplifies as follows: 

dp{t,T) ^H{r) \{rt + htLt)v{t,T)dt + h{t)d(l){t) \ + v{t' ,T)dH{t) 

= {rt + htLt)p{t- ,T)dt + H{t-)h{t)d(j}{t) - v{t~,T)dH{t). (57) 

Let us now try to find the dynamics of 0(t). Since — ^' {%l){t),Xt) , 
where V'(i) is given as in (19), Ito's formula leads to 

dm = -(rt + /i,Lt)e-^*('-"+""^")'^" {m,Xt) dt 



{^,xt) + {mAA%x,)'^ dt 
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where we had used the difFerentiabihty proved in Lemma A. 2 and the semi 



martingale representation formula of our Markov chain Xt given in Eq. (13). 



As (f>{t) is a (F, Q)-(local) martingale, its drift term is zero and, therefore, we 
obtain the dynamics 

d^{t) = e-/o('-"+''"^")rf«(V7(t),dMQ(t)) (58) 

and all together, we get 

dp{t,T) = {rt + htLt)pit~,T)dt + H{t-){^{t),dM^{t)) - v{t- ,T)dH{t). (59) 

we have dH{t) ^d^f + {l~ 



A.l 



By the Doob- Meyer decomposition of Lemma . 
H{t~))htdt. The last step follows from the fact that on the event t > t, we have 
v{t-,T) = p{t-,T), by definition of p(i,T), and also p{t- ,T) = (V-W,^*-)- 



Therefore, we can write the pre-default risk-neutral dynamics of the bond as 
dp{t,T) 



in + h,{L,-l))dt + ^^^^^^^-dif 



(60) 

□ 



Proof of Proposition \3.S\ The bond price dynamics under the real- world mea- 



sure follows directly from plugging Eq. (|14|) into (|60|) to get: 
dp{t,T) 



rt + ht[Lt - 1) + ^ , , , I dt 



{m,dM^{t)) 



(20) since 



where we also used that ^ = (see Lemma A.l ). Finally, we get the dynamics 

(61) 



{m. {A'{t) - {A^y{t))X^) {{A{t) - A^{t))m,X^) 



When Xt — e^, the previous quantity equals 



A(t) :=EK..W-«?,W)|^. 



and (61 ) can be written as ({Di(t), . . . , Z?Ar(t))', Xt) 



□ 



B Derivation of the generator of {t, Vt, Ct, Ht) 

We start by changing our notation to be more consistent with the framework in 



Bielecki and Rutkowski (2001). To this end, let 

AT 

Ct :=Ell{x,=e,}- 



(62) 



27 



Note that {Ct)t is a Markov process with values in {1,...,N} and infinites- 
imal generator A{t) = [ai.j(t)]i.j=i,..._Ar. In particular, for any function g : 
{0,...,7V}^M, 

MS{t) giCt)-g{Co) - f\Ag){C^,u)du, 



Jo 

is a martingale under P, where we have used the notation 

N 

^5(^,0 :=[g'^'W]» = E«^^ W50')' ^ith g - (.9(1), . . . , g(iV))'; 

3 = 1 



c.f. Proposition 11.2.2 in Bielecki and Rutkowski (2001^. Note that this result 



follows directly from the semimartingale decomposition (111) by multiplying (from 
the left) both sides there by g'. In particular, also note that 



MS{t) = g'Ar{t), and Mj{t)^Hl^ a^.^^{u)du, (63) 



where we used notation ( 11 ). Let us assume that V admits the following Markov- 
modulated dynamics: 



N 



dVt = a^^ dt + dWt + Pc^_ , dMjit) - 7o^. dft 



(64) 



where ai{-, •, z), i3i(-, •, z), •, z), "fi{-, •, z) are deterministic smooth functions 

in [0,00) X M for any i,j € {1, . . • , N} and z G {0, 1}, and all the coefficients 
in (64) are evaluated at {t,Vt- , H{t^)). In terms of the processes (lO)-(lll and 
using (63), we first note that 



N 



N 



(65) 



N 



N 



= ^ P.,,it,Vt-,H{t-))^c,^=^}dHi- J2 ^.j{t^V^-M{t-))Hla^^^^{t)dt 

N N 

= E E(^^^^- - ^^^)(^' ' H{t-))dHl ' - J2 Pc,., (t^ Vt- , (t-))a^„, {t)dt. 

i=l j=^i j=l 

and, in particular, 

N 

= E E ^t" ' H{t-))AHl'^ - 7c^_ {t, 14- , H{t-))AH{t), (66) 

where = ftj Next, let f{-,-,i,z) e C^'^{[0,oo) x M), for each i = 

1, . . . ,N and z G {0, 1}. We want to find the semimartingale decomposition of 
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f{t,Vt,Ct,H{t)). Applying the Ito's formula (seeing Ct as simply a bounded 
variation process), we have that 

fit,Vt,Ct,Hit)) = f{0,Vo,Co,H{0))+ f ft{u,Vu,Cu,H{u))du (67) 

N 



1 



'dcJvdWu + - / hv^lju 



J2 {/KK,C„,iJ(u))-/(u,K-,C„-,i7(u-))} 



0<u<t 



Since r is not a transition time of C a.s., we can write the last term in the above 
equation as follows: 



Jt-= {/KK,C„,o)-/(«,K-,C„-,o)} 

0<u<t/\r 

+ {fir, Vr- - 7c^_ (r, K- , 0), C'r- , 1) - f{u, K- , ^r- , 0)} 
+ ^ {/(w,K,C„,l)-/KK-,C„-,l)} 



t/\T<U<t 



+ ^ {/(", K- - 7c„_ K- , 0), , 1) - f{u, , , 0)} di7(u). 
Next, using the local martingales ^ and (10), we have 

Jt = Y.Y. / [/KK-+/3°,,j-,i?(^-))-/KK-,*,i?(w-))]rfAf;'^' 
+ / -7,^_(u,K-,o),c„-,i)-/KK-,c„-,o)}dcj: 



+ / E /(«,K + /3;! ^,j,i7(u))-/(u,K-,C„,i7(u)) du 

+ / {/(w,K-7c„(":K,0),C„-,l)-/(u,K,C„,0)}(l-i?(M))/^^^d?/, 
Jo 

where we had also used that = V^- , H{u) — H{u~), and C„ = C„- a.e. 
and, hence, the integrands in the last two integrals with respect to du can be 
evaluated at {Vu,Cu,H{u) instead of , C„- , iJ(w~). All together, we have 
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the semimartingale decomposition 



fit,Vt,CuH{t)) = fiO,Vo,Co,HiO))+ Cf{u,Vu,Cu,H{u))du + Mt, (68) 







where {Ait)t is a local martingale and Cf{t,v,i, z) is the so-called generator of 
{t,VuCuH{t)) defined by 



dt dv 



N 

ai{t,v,z) - ^I3i^j{t,v,z)a,^j{t) + {1 - z)hij, 
^\ pa 



2 dv^ 

+ ^o.ij{t) {f{t,v + f3lj{t,v,z),j,z) - f{t,v,i,z)) 

+ {fit, V - 7,(t, V, 0), i, 1) - fit, V, 1, 0)} (1 - z)h,, (69) 



for each i = 1, . . . , N . The local martingale component in (68) takes the form: 

N 



+ / {/(«,K.- -7.(^^,K-,0),i,i)-/(w,K-,«,0)}i{c„_=z}< 



+ i ^9.^(w,K,«,i?(u))i?;dWn|, (70) 
where the functions l3i,j, j3'^,, and di are evaluated at iu,Vy^~ , H{u^)) and we 



used the notation (111 



C Proof of the verification theorems 



Proof of Theorem 5.1 We first note that in the post-default case, the process 



( 26 ) takes the form 



dV^''-'^ = V:-'^^ I [r^^ + (Mc. - ^c. )] ds + TTsa^^ dWs | , (71) 
Ft"'*'" =v, {t<s< R). 
Define the process 

M; := w{s, a), {t<s< R), (72) 

for an admissible feedback control tt^ :— tt^._^{s,VJ''^''") £ Ativ,i,l). For sim- 
plicity, through this part we sometimes write Vu or instead of 1/^'*''". We 
prove the result through the following steps: 



30 



(1) By the semimartingale decomposition (68), it follows that 

M; = + ^ Riu, V:, C„, Tr^)du + Ms ~ Mu 

where 



Ms 



N . ^s 



R{u,v,i,TT) = w^{u,v,i)+w^{u,v,i)v{ri+Tr{fii~ri)) 

I 1 2 2 2 

+ E"«j(") ^Ml{u,v,i) -w{u,v,i)) . 



(73) 



(74) 



We have that R{u, v, i, tt) is a concave function in tt since, by assumption, w^^ < 
0. If we maximize R{u, v, i, tt) as a function of tt for each (u, v, i), we find that 
the optimum is given by (40 1. This imphes that 



R{u, V, i, tt) < R{u, V, i, 'Ki{u, v)) = w^{u, v, i) + rivw^iu, u, i) — rji 



wl{u,v,i) 



where the last equality follows from Eq. ( |37| ). Next, let us introduce the stopping 
times Ta^b inf{s > t : Vs > b^^, or Vg < a}, for fixed < a < v < < oo. 
Then, using the notation Et[-] = E[-|tJt], we get the inequality 



N 



<M^+ E 



{w{u,V^-,3)-w{u,V^-,i)}dM:;' 



N 



■EE. 



with equality if tt = tt. Since 

sup |w(u,K,«)| < ^1, sup \VuT:i{u,Vu)w^{u,Vu,i)\'^ < B2, 



t<u<ra faA-R 



for some constants -Bi, i?2 < 00, we conclude that 
with equality if tt = tt. 
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(2) In this step, we show that 



Urn Et 



where = ^{s, F/'*'", C^). Note that (|38[i) imphes 

2 



E 



< Bi + B2E 



(75) 



for some constants -Bi,i?2 < co- Next, we note that tt satisfies (39) since 



\VTT,{S,V)\ 



<G{s){l + v), 



in hght of (38-ii). Hence, we can apply Lemma C.l below (with tt^ = 0) and 
obtain 



sup Et 

0<a<ii<fc-i<oo 



2Et 



sup (vf~Vl 

,t<s<R ^ 



< 00. 



Using Corollary 7.1.5 in Chow and Teicher (1978), we conclude (75) 



(3) Finally, if w_ is non-negative, then Fatou's Lemma implies that 

Et [U{V^)] = Et 

< hm inf Et w{R A t^.^, V]^^ , Cr;,^^ , ) 

a.o— :-0 L ■ u.o , 

<w{t,v,Ct)^Et \u{Vi) 



lim inf w{R A Ta,6, V^j^ , Cj?^Ar„ i, ) 



for every admissible feedback control tt^ = tt^ (s,V7^'*'") G ^f(w,i,l). For a 
general function w (not necessarily non- negative), we proceed along the lines of 
step (2) above to show 



lim Et 

a,b->-0 



for any feedback control tt^ — tt^ (s, V/'*^") € ^t(w,*, 1) satisfying (39). □ 



Proof of Theorem \5.S^ We prove the result through the following steps: 
(1) Define the process 

:= w{s, - H{s)) + w{s, V^^'^\Cs)H{s), (76) 



where V^' is the solution of Eq. (26) for an admissible feedback control 



TT, := (7rf,^f) (^^ (s, Fl'*'"), < (s, V;j*''')) € A(^;,^,0). 
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For simplicity, we only write = F/'*'". Using the same arguments as in 
Eq. ( 65 1 and the decomposition Q , the process ( 26 1 can be written as 



ds 



•TT^ a^dWs + irP {l~H{s-))dH{s) 



(77) 



N 



(l-ff(,s-))^^vr. 



for s G {t,R), with the initial condition y^'^^*''" = where 9i is defined in (34). 
By the semimartingale decomposition ( 68 1 with the coefficients given by ( 32 1 
and 



f{s,v,i,z) := w{s,v,i){l — z) + w(s, v,i)z, 



it follows that 



M; = M[ + ^ R{u, IC, Cu, irl^ , , H{u))du + Ms - Mt 



Here, A^g — is given by 
A' . 

eIe 

i=i ^ 



ipj{u) - tpi{u) 



-P Vj 



i>i{u) 



j - w{u,Vu-,i) 



dMl 



where fff {u, v) (1 — H{u))Tif {u, v) and 

w{u, V, i) — w{u, V, — H{u^)) + w{u, v, i)H{u^). 



(78) 



Similarly, i?(w, u, i, tt'^, tt-'^, 1) is defined as in (74 1 with tt = tt"^, while 



R{u,v,i,'K ,Tr ,0) := Wu{u,v,i) + Wv{u,v,i)v\ Ti + t: (fii — ri) + n 9i{u) 



3¥=i 



W [ U,V [ f + TT 



p yjj 



hi [w (u, W (l — TT"^) ,ij — w (m, V, i)] 



,j ] - w{u,v,i) 
(79) 



Note that, under our assumptions, {tt^jTt^) — > R{u,v,i,TT^ ,1) admits a 
unique maximal point (tt"^, tt^) for each {u, v, i), since (i) i?^s^s < 0, (ii) R^^p^p < 
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0, and (iii) Rt^s^p = 0. Indeed, (i) follows from our assumption that w„„ < 0, 
while (ii) is evident from the calculation 



V Wyy I U,V { 1 + TT 



p 



+ hiv'^W.^^ (u, v[\- TT^) , i) 



and the fact that w^^ < 0. The optimum is given by Eq. (44 1 with p defined 
imphcitly by Eq. (liTl). In light of Eqs. (l37|-(l^, 



Let 



R{u, V, i, TT^ , ,z) < R(u, V, i, 7rf (u, w), 7rf (m, v), z) — 0. 
Ta,b ■■='vii { s>t ■.Vs> or Vs < a, or 7r^_ > 1 - a, 



or 7r„ 



< a — 1 , for some j ^ i 



for a small enough a,b > 0. Using similar arguments to those in the proof of 
Theorem 1 5.1 1 we can show that 

< = wit, V, Ct)(l - H{t)) + wit, V, Ct)Hit), 



with equality if tt'^ — and tt^ — . 
(2) We now show that 



lim Et 


wiR A Ta.b, Vl^,^^^,CR^r^,)HiR A r,,6) 


= Et 


C/(y|)ff(i?) 


(80) 


lim Et 

a.b-^a 


A Ta.fc, V/^Ar„,, . CRr,r^ ^)HiR A ra,fc) 


= 


UiV^)HiR) 


, (81) 



where 7J(s) = 1 - i7(s). Note that ( [SOpT] ) will imply that 



wit,v,Ct)Hit) + wit,v,Ct)Hit) ^ lim Ef 

a,6— >0 



Ef 



We only prove (80) ((81 1 can be treated similarly). Note that (38-i) implies 



E 



wiRATa,b,V^^^^ ^,CRAT^ ,)HiRATaM) 



<Bi + B2E 



for some constants -Bi, i?2 < oo- Next, we note that satisfies (391 since both 
w) and w satisfies (38 -ii) by assumption. Also, satisfies (28) by assumption. 
Hence, we can applying Lemma |C . 1 1 below and obtain 



sup Et 

0<a<ii<fc-i<oo 



V, 



< 2 



+ 2Et 



sup 

,t<s<R 



< 00. 



Using Corollary 7.1.5 in Chow and Teicher (1978), we conclude (80). 
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(3) Finally, if w is non-negative, then Fatou's Lemma implies that 



E, [U{V^)] - E, 



liminf w{R A Ta,6, V^^ .,CRATa b)H{R A Ta^b) 



- Et liminf A Ta,b, VSat^^^Cr^t^ ,)H{R A Ta,b) 



< liminf Et 

a,b-^Q 



W{R A Ta,b, V^j,^^^^,CR/,r^ jH{R A Ta,b) 

+ w(R A Ta^b, V"flAr„ ,,Cbat^.JH{R A r^^h) 



< w(t, Ct)H{t) + w{t, V, Ct)H{t) = Et U{V]^) 



for every admissible feedback control tIs — i^c _ is,V^^*'" , H{s )) € ^t(w,'i,0). 
For a general function w (not necessarily non- negative) , we proceed along the 
lines of step (2) above to show 



Hm Et 

a,b-^0 



j{R A Ta,b, ^flAr„ . , C'flAr^, Ji? (i? A Ta.fc) 



+ wiR A Ta,b, V^flAr. . , CrAt,,,)H{R A T^.h) 



= Et [[/(V-;?)] 



for any t-admissible feedback controls tTs = tt^ (s, VJL'*'", )) e ^t(^^,*,0) 
satisfying ^ and (p8|. ° □ 



Lemma C.l. Let 7r,f (s, w, z) and irf (s, u, z) be functions such that (H) admits 
a unique nonnegative solution (V7^)sg[t.i?j • We also assume that tt^ satisfies 



39) and tt^ satisfies (28). Then, the solution of (ID satisfies the moment 
condition: 



Et 



sup \v:~vr^ 

t<u<R 



(82) 



for some constants Ci , C2 < c» . 

Proof. For simplicity, we write Vs instead of . Let us start by recalling that 
we can write ( 77 1 in the form 



N 



dVs = + ^c, + 51 '^c _ (s) - 7c _ '^Cf , (83) 



taking the coefficients as in (32). Due to (28) and (39), we can see that 
|ai(s, u, z)\ < E{s){l + v) for a locally bounded function E. Hence, by Jensen's 
inequality and the previous linear growth. 



{u,Vu,H(u ))du 



< KT r + rVt + 



VA^du 
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for any s G [t,R], where t := R — t and k denote a generic constant that may 
change from hne to hne. Shmlarly, denoting = inf{s > t : \Vs\ >h}{h> v) 
and using Burkholder-Davis-Gundy inequahty (see Theorem 3.28 in [Karatzas 



and Shreve| ( |1998| ) or Theorem IV. 48 in |Protter| ( |1990[ )) and Jensen's inequahty, 



Et sup 

t<S<Rf\Tb 



{u,Vu,H{u))dWu 



< ktE / |crp^7r„Ki| du. 



We can then again use ( 39 1 to show that 



Et sup 

t<s<RATt, 



RATb 



< KT{T + T\VtY+ / \Vu-Vt\'du). 



Next, using again Burkholder-Davis-Gundy inequahty (see, e.g.. Theorem 23.12 
in 



Kallenberg (1997)), 



Et sup 

t<S<RATb 



-f^ _iu,V,-,H{u-))de, 



< KEt 



RATb 



V^-7T^\^dH{u). 



Using (281, 



E, 



/RATb pRATb 
\V^-TrPfdHiu) < K,Et \V^-\''dH{u 



KEt 



RATb 



\Vu\''K^{l - H{u))du. 



Then, we can proceed as before to conclude that 



Et sup 

t<S<RATb 



7c _{u,V^-,H{u-))dil 



j-RATb 

< kt{t + T\Vt\'' + \Vu'Vt\'du) 



Using a similar argument, we can also obtain that 

2 



Et sup 

t<S<RATb 



2 I- RATb 

f3^^^_^iu,V^^,H{u-))dMjiu) <KT{T + T\Vt\^ + J^ \Vu-Vt\ 



Putting together the previous estimates, we conclude that the function 7f,(r) := 
Et supt<s<r/\T-j^ \Vs — VtP can be bounded as follows: 7&(r) < kt{1 + v^) + 
K 'yb{u)du. By Gronwall inequality, we have 



76 



(i?) < A«(i?-t)(l + w2)e"(^'~*), 



and (82 1 is obtained by making 6 — oo. 



□ 
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D Proof of Explicit Constructions 



Proof of Lemma 6.2 For fixed i and s, consider the function 



f{pi,s,i) := 0i{s) 



We first observe tfiat f{pi,i,s) is a continuous function of pi in the interval 
(Mj, 1). Indeed, we can write the above summation as 



E 



E 



and since 1 < when < we have p, + < for 

Pi S (Mi, 1). Moreover, the previous decomposition also shows for each fixed s, 
f{Pi,s,i) is strictly decreasing in pi from (Mj,l) onto (— cxi,oo). This implies 
the existence of a unique Pi{s) such that fipi{s), s, i) — 0, for any s > 0. In light 
of Kumagai (1980) implicit theorem, we will also have that Pi{s) is continuous 
if we prove that f{pi,s,i) is continuous in (j)i,s). The latter property follows 
because, by assumption, Uij and afj are continuous, which implies directly the 
continuity of the functions 9i. The continuity of the functions tpj will follow 
from a similar argument to that of Lemma |A.2[ □ 



Proof of Proposition \6.3\ 
(1) It can be checked that the function (p^{t,v,i) 



Dirichlet problem (37l, if the functions K{t,i), i 



- log(w) + -ftr(t, i) solves the 
1, . . . , satisfy the system 



given by Eq. (46), which may be written in matrix- vector form as 



(84) 



As 



is continuous in 



where [F{t)]i j = —aij{t) and b{t) = — 
[0, T] by hypoth esis, we have that the system admits a unique solution in 
[0, R] by Lemma 



6.1 



where R < T. Moreover, ip^{t,v,i) G C°2 due to con- 
cavity and increasmgness of the logarithmic function, and, under the choice 
D{t) = maxi^i ..._jv i4r(t, z) and G{t) = 1, the function (p^{t,v,i) satisfies the 
conditions in (38). Therefore, applying Theorem 5.1 we can conclude that, for 
. , A^, 'P^it, V, i) is the optimal post-default value function. 



each i = 1, 



(2) Plugging the expression for ip^{t,v,i) inside Eq. (40), we can conclude im- 



mediately that fr^ {t, j) 

(3) It can be checked that the vector of functions [Tpf'{t,v, 1), . . . ,Tpj^{t,v, N)], 
where ^f {t,v,i) — log(w) -I- J{t,i), and the vector p{t) solving the nonlinear 
system of equations ( 45 ) simultaneously satisfy the system composed of Eq. (|4T|) 
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and Eq. ^ if the vector of functions J(t) = {J(t, 1), J(i, 2), . 
the system (|47|), which may be written in matrix form as 



.,J{t,N)) solves 



Jt{t) = F{t)J{t) + b{t), J{R) = 



-Ci^At) (J^^)^ iFit)]'^ 



hit) = - r, - Mt)9,{t) - h, (log(l - p,{t)) + K{t, i)) 



^a,j(i) log(l+p,(t) 



It remains to show that such solution vector J{t) is unique. As we know that 
Pi{t) is continuous in t, then we have that each entry [F{t)\i_j is a continuous 
function of t. Moreover, the vector K{t) is continuous, as it solves the system of 



differential equations given by (84). As af^ is continuous by assumption, then 



'i/'j's are continuous, thus 9i is continuous, and consequently hi{t) is continuous 



in [0,i?]. Using Lemma 6.1 we obtain that the solution vector Jit) must be 
unique. Moreover, under the choice of D{t) — maxi^i... jv J{t, i) and G{t) = 1, 



we have that ip^(t,v,i 



satisfies the conditions in (38) 



As the logarithmic 

function is increasing and concave in then (/jf (t, v, i) G Cq'^, therefore it must 
be the optimal pre-default value function by Theorem |5.2| □ 
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